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Abstract 

Computing the trajectories of particles in generic beamlines is an important ingredient of experimental 
particle physics, in particular regarding near-beam detectors. A new tool, Hector, has been built for such 
calculations, using the transfer matrix approach and energy corrections. The limiting aperture effects are 
also taken into account. As an illustration, the tool was used to simulate the LHC beamlines, in particular 
around the high luminosity interaction points (IPs), and validated with results of the Mad-X simulator. 
The LHC beam profiles, trajectories and beta functions are presented. Assuming certain forward proton 
detector scenarios around the IPS, acceptance plots, irradiation doses and chromaticity grids are produced. 
Furthermore, the reconstruction of proton kinematic variables at the IP (energy and angle) is studied as 
well as the impact of the misalignment of bcamlinc elements. 
PACS: 29.27.-a ; 41.85.Ja 
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1 Introduction 

In the context of the LHC beamlines and in particular of very forward, near beam detectors, we introduce 
a new fast simulator, called^ Hector, for the transport of single particles through generic beamlines, 
although primarily dedicated to the LHC. The simulator is based on a linear approach of the beamline 
optics, implementing transport matrices from the optical element magnetic effective length, and with 
correction factors on magnetic strength for particles with non nominal energy. Hector deals with the 
computation of the position and angle of beam particles, and the limiting aperture of the optical elements. 
It has been designed to be fast, light and user friendly and its object oriented structure, using the ROOT 
framework [1], helps its usage, its maintenance and the future improvements. 

In forward physics, the measurement of the position and the angle of particles in dedicated detectors 
(as so-called roman pots for example) hundreds of meters away from the interaction point (IP), helps in 
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reconstructing the kinematics of the event at the related central detector. Hector links the information 
from very forward detectors (VFDs) to the one from the central detector, by precise calculation of the 
particle trajectory. 

We present here the physics content and the simulation techniques used in Hector, as well as its 
overall performance. The validation of Hector output comes from direct comparisons with the Mad-X 
simulator [2], for the beta functions and relative positions of both LHC beams. Some physics results, like 
beam profiles and irradiation levels due to the forward proton from inelastic pp pX interactions are 
plotted for two reasonable scenarios for the very forward detectors. Finally, first studies of the forward- 
scattered proton kinematics reconstruction are performed, resulting in the discussion of limiting factors 
and best strategies for the VFD designs. 

2 The LHC beamlines 

Even if Hector's class structure is compatible with the most frequent beamline elements, the physics 
motivation on backstage is driven by the forward physics at the LHC. This includes detection of protons 
which have for example been quasi-elastically scattered by exchange of a photon. The reconstruction of 
the proton momentum using information from a VFD requires a good knowledge of the LHC beams and 
related beamlines. Such a beamline is described here below ; note that both LHC beams are basically 
identical. 

2.1 Beamline elements and beam parameters 

The position and lateral size of the LHC beams are mainly controlled by two types of elements : 

• focusing (quadrupoles) or deflecting (dipoles, kickers) magnets; 

• collimators and optical element apertures, limiting the lateral beam extension. 

Quadrupole magnets can focus the beam either vertically or horizontally. Dipole magnets also come 
in two kinds : rectangular dipoles with a straight shape, and sector dipoles which are bent to match the 
beam curved trajectory. The former are mainly used in the long and short straight sections (LSS, SSS) 
of the beamline, while the latter are placed in the dispersion suppressors (DSs) and the bending sections 
(ARCs) [3], [4]. The kickers are dipole magnets dedicated to produce the crossing angle at the IP. 

Some parameters are of direct interest when describing the nominal 7 TeV proton beams: 

• the horizontal and vertical emittances ex^y ^ 500 x 10^® /im, which are constant along the orbit 
according to Liouville's theorem ; 

• the lateral beam size at the IP cr = 16.6 fim ~ \/e/3* , where j3* is value (= 0.55 m) at the IP of the 
so-called beta-function ; 

• the lateral vertex size (or the transverse size of collision region) a" = 11.76 /im = a j \f2 ; 

• the angular divergence at the IP ue = 30.2 /irad = y^e//3* ; 

• the crossing angle is equal to 2 x 142.5 /irad, in vertical (horizontal) plane at the IPl (IP5) ; 

• the energy spread ue = 1.129 x IQ^'^E. 
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2.2 Beamline map 



From here onwards, focus will be put on interaction region 5 (IRS). The elements affecting the beam will 
be described starting at interaction point 5 (IP5), going eastward when looking from the center. This is the 
propagation direction of beam 1. The other side of the beamline is very similar. As the region of interest 
for the forward physics related detectors extends up to approximately 500 m from the IP, our description 
of the beamline stops there. Within this part of the beamline, all aperture shapes are rectellipses [4]. 

Before 85 m from IPS 

After the IP, beam 1 gets separated from the incoming beam 2 because of the crossing angle. The two 
beam paths stay close (~ 4 mm) for the first 85 m. They encounter the first three collimators at 19 m, 
45 m and 55 m from the IP (the so-called low-/? inner triplet), then a kicker (30 m) and then quadrupoles 
(H: 23 and 47 m ; V: 32 and 39 m) and consecutive dipoles (R-Bend) from 60 to 81 m, which deviate the 
two beam paths to lead them to their final separation of 194 mm. The dipoles are positioned symmetrically 
with respect to the two beams. 

From 85 m to 269 m from IPS 

At about 160 m from the IP big R-Bend dipoles are used to make the beam orbits parallel, so from there 
on the beams get to their nominal relative positions. These dipoles are positioned in such a way that the 
beam 1 (2) enters (exits) it perpendicularly. Four more collimators are located respectively at 141, 149, 
184 and 256 m from the IP5. Two kickers, at 165 and 199 m, lead them to the ideal path that they would 
have followed if there were not any crossing angle, due to the so-called closed bump orbit. The beam is 
kept under control with five more quadrupoles (H: 194, 260, 264 ; V: 168 and 226 m ) and one dipole 
(153 m). This is the end of the straight section. 

From 269 m onwards 

There the beam path starts its bending. The beamline is composed in alternations of quadrupoles (> 6) 
to keep the beam focused and sector dipoles (> 8 S-dipoles) to bend it. The first sextupoles with non-zero 
field occur only after 439 m from the IP. 



3 Simulation techniques 
3.1 Physical description 

The simulator relies on a linear approach to single particle propagation. A beamline consists of a set 
of optical elements, amongst dipoles, quadrupoles, drifts, collimators, kickers and VFDs. The optical 
elements are described by their magnetic field, their length and their aperture. In turn, a set of particles, 
with all possible smearings of initial positions, angles or energies, is propagated by Hector through the 
beamline, particle by particle. In electromagnetism, the influence of an external magnetic field is given 
by the Lorentz force. Let's consider the Taylor expansion of the vertical component of magnetic field By, 
around its central value : 

e ^ , , e ^ e dBy 1 e d^By , 

-By{x) = -By + -—^X + --—^X^ + . . . 

p p p OX 2p ox^ 

where p is the momentum of the particle and e its electric charge. The terms of this sum are interpreted 
as respectively dipolar (-^ = p^y)^ quadrupolar (fc = ^-g^) and sextupolar fields. In the co- moving 
coordinate system, neglecting small deviations (x R, y ^ R) and small momentum loss ( Ap <C p) , this 
leads to the following equation of motions for a particle traveling along path length s through a magnetic 
element [5] : 
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\ y"{s) + k{s)y{s) = 0. 

The solution (a;(s), a;'(s), j/'(s)) to these equations can be expressed as a hnear combination of 
the initial values (zq, a^Q, j/Oi J/o)' where the rotation matrices are defined by the properties of the optical 
element (length and magnetic field). The typical values R ~ 200 m and x < 0.01 m match the first 
assumption (x <C i?). As Ap <C p is not always verified, we apply a further correction to the magnetic 
field (for mode details, see section 3.2) of each optical element of the beamline. 

In other words, for the simulation of the transport of a particle in a beamline, each beam particle is 
represented by a phase space vector and each optical element by a transfer matrix by which the vector is 
multiplied. The propagation of a single particle is thus the rotation of the phase space vector by the n 
optical element matrices. 



X{s) = X{0)MiM2...M^ 

The whole beamline is modeled as a single transport matrix acting on each particle phase space vector 
(no intrabeam interactions). The optical element description also refers to its physical aperture. When 
a particle is propagated through an optical element, two tests check whether its path is compatible with 
the element acceptance or not, at its entrance and exit. The fact that a bent path could hit a central part 
of the element is neglected. A particle which does not pass through the element acceptance is flagged, 
allowing to know the stopping element for given particle parameters (for the computation of acceptance 
and irradiation rates for instance). 

3.2 Implementation 

The simulator has an object-oriented design, using the C++ ROOT framework [1], with dedicated classes 
describing : 

• the beam particles. The 6-components phase space vector is X = (x, x' , y, y' , 1), with horizontal 
{x,x') and vertical {y,y') coordinates and angles ; E is the particle energy. The sixth component is 
just a factor used to add an angular kick on the particle momentum direction. 

• the optical elements (dipoles, drifts...) inheriting from a common class. 
The 6x6 transport matrices can be decomposed into blocks 

/ A A 00\ 

^ ^ 

B B 

M„nits- Q B B 

V V 10 
\ K K K K 1 J 

where 

• A (and B) blocks refer to the action (focusing, defocusing, drift) on horizontal (and vertical, rcsp.) 
coordinate and angle. 

• V terms reflect the dispersion effects of the horizontal dipoles on off-momentum particles. 

• K factors arc the angular action of kickers. 
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A beamlinc is the implementation of a list of optical elements with such transport matrices. The 
introduction of dispersion terms allow a proper description of off-momentum particles {Ap/p ^ 0). The 
dispersion function can be defined from eq. (1), for horizontal dipoles (fc = 0), taking Ap/p ~ 1 : 

+ ^OM 4^ (2) 

Developing the solution of this equation leads to a correction term for the deflection of off-momentum 
particles in the dipoles : 

Ap 

2^off-momontum(s) = x{s) + D{s) . (3) 

P 

For explicit description of these V terms, refer to the appendix. Eq. (9). Besides the dispersion 
correction which is valid for low Ap, the actual dependence of the particle deflection on its momentum 
is also taken into account. This is the chromaticity, or the energy dependence of the transport matrix, 
implemented by rescaling the magnetic field terms {R, k, K) with a factor ( J'^^ ) ■ The propagation of 
particles different from protons is also possible by rescaling these magnetic field terms : 

ki{AE,mp,qp) = — — , h = R,k,K; (4) 

PO - Ap ^proton 

where qp and rup are the particle charge and mass, respectively. The ROOT framework easily interfaces 
Hector with the output of the high energy physics event generator Pythia [8] . Forward particles from 
the final state can be then propagated through the beamlinc via Hector, for example in photon induced 
interactions. 

As Hector has been designed for the LHC bcamlincs, one class can parse the official LHC optics 
tables, interfacing them directly into the code at run-time. For a better compatibility with different table 
layout. Hector can recognize column headers and is thus compliant with different table types [9]. All the 
default parameters are hard coded (in H_Parameters .h), like nominal beam energy, position and angle, 
and their divergences, as well as the nominal crossing angle, but these default value can obviously be 
rewritten at run time. The effects of the kickers or the bending of the sector magnets can also be easily 
switched on or off. For the following analysis, LHC optics version 6.500 has been used, corresponding to 
/3* = 0.55 m at the IPl and IP5. 

Limitations : sextupoles and magnets of higher order (multipolcs) are neglected in Hector (harmless 
the first 430 m after the IP). The energy dependence of the transfer matrices is taken into account, but 
neither beam-beam interactions, nor the field non-uniformities. The beam particle mass is also neglected 
compared to its energy. 



4 Validation and performances 

4.1 Performances 

Hector has been designed to be fast and light. All the code and its library altogether weight less than 
2.5 MB for about 12000 lines of code. The figure 1 depicts its performance, leading for a first estimate to 
a CPU time of 4 /is per beam particle per beamline optical element at the LHC. As expected, the CPU 
time is linear with the number of optical elements and the number of protons. 

4.2 Validating Hector with Mad-X 

In order to validate Hector, cross-checks have been made with the output of Mad-X [2], such as the 
beam beta functions (Fig. 2) and the beam relative positions (Fig. 3). For each point, the beta functions 
are computed from the beam omittance and the RMS of its profile in the transverse plane. The beam beta 
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Figure 1: Hector speed performance : CPU time with respect to the propagation distance {left) and the 
number of optical elements {right) in the beamline, computed for a 10^ particle beam around the IR5. 
As the propagation of particles corresponds to the multiplication of matrices, the computation time is 
directly proportional to the amount of optical elements in the beamline, as well as the number of beam 
particles. For a time estimation for large number of particles and long bcamlincs, a good approximation 
is 4 us per optical element (including drifts) per particle. For this test (3.43 /US per element per particle). 
Hector is run on a regular laptop with 1.7 GHz Intel Pentium Ccntrino processor. Linux rel. Ubuntu. 



functions, which describe the lateral beam size variation along the orbit, as output by Hector, match 
perfectly the results from Mad-X. The lateral beam position shown in Fig. 3 is relative to its ideal path, 
defined for disabled kickers (and no crossing angle). These graphs (Fig. 2 and 3), and their equivalents 
in backward direction, arc the strongest tests for Hector, where smearings of position, angle and energy 
are applied as previously explained. 

4.3 Validating off-beam energy simulation 

The comparisons of Hector and Mad-X results directly validate only the simulation of a 7 TeV proton 
propagation, and additional checks are necessary to validate simulations of the off-beam energy particles. 
This has been done by comparing numerically the exact calculations and the Hector predictions for a 
couple of proton trajectories at energies between 10 and 100% of the beam energy. It was found that for 
all the energies of interest for the VFDs, i.e. for energies above 80% of the beam energy. Hector remains 
very accurate, providing the particle lateral positions with precision at least at a few micron level. 
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P functions - beam 1, forward 



P - Hector 
- Hector 
P^ - MAD-X 
p' - MAD-X 




P functions - beam 2, forward 



P - Hector 
p' - Hector 
P^ - MAD-X 
P' - MAD-X 




Figure 2: Beta functions (3^ and Py for beams 1 (above) and 2 (below) near the IP5, going forward: 
these functions depict the beam size variation, and are computed from the beam emittance and the beam 
lateral profiles at successive positions. Indexes x and y are corresponding to respectively the horizontal 
and vertical directions. The plain and dashed curves, coming from Mad-X, match exactly Hector's 
output (squares). 
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Figure 3: Relative horizontal positions for beams 1 (above) and 2 (below), around the IPS, in the forward 
direction. The relative position is defined as the difference along the ideal orbit between the beam lateral 
positions with and without the crossing angle. This relative position graph (squares) matches Mad-X 
predictions (plain lines). One clearly sees the crossing angle at the IPS (s = m) and the effects of the 
kickers that move the beams back to the ideal orbit. 
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5 Beamline simulation 



By using the HECTOR classes and library physics studies can be performed by means of the particle 
propagation inside the beamlines, away from the interaction point. 

5.1 Trajectories 

Knowing the optics tables for both LHC beams, their trajectories can be compared simultaneously, in 
both top and side views, for the two LHC beams aside (incoming beam 2 and outgoing beam 1, at the IPS 
and IPl: Fig. 4). The top view shows the beams on the horizontal plane, clearly depicting the crossing 
angle at the IPS, and the beam separation after 70 m away from the interaction point. The bending of 
the sector dipoles has been switched off in order to makes graphics more clear - this is why both beams 
are straight and parallel after 2S0 m. However, the optical elements have been shifted (without tilt) in 
the horizontal plane by the half of the beam separation distance, from 180 m away from the IP, in order 
to match the ideal beam path: a proton with nominal energy and on the ideal orbit should travel through 
optical elements in their geometrical center. The side view in turn shows the beams in the vertical plane, 
emphasizing the difference between IPl and IPS. In addition, the major optical elements have been drawn: 
rectangular dipoles in red, sector dipoles in light green, and quadrupoles in yellow and blue. 

5.2 Beam profiles 

Lateral beam profiles along the beamline are of interest, for instance for studies of very forward detectors 
(Fig. 5). The beam size and shape are obtained after propagation of beam particles using initial dispersion 
in position, angle and energy at the interaction point. Assuming a Gaussian distribution for these variables, 
the evolution of both beams can be seen in the transverse plane for example at 220 m and 420 m away from 
the IP. As expected, the beam particles are distributed symmetrically on the vertical plane (y = 0mm), 
and are already separated in the horizontal plane {x = ±97mm) . The symmetry between beam 1 and 2 is 
striking. The 3cr contours (in red) guides the eye for the lateral beam extension. In the angular (9x,0y) 
plane, the beams evolve according to the focusings and the defocusings they undergo, even if the mean 
value of these angles remains effectively equal to /irad. The evolution of the beams in the phase space 
has also been drawn (Fig. 6). Once again, the similitude between both beams is clear, even if they differ 
from each other by slight details. All these profiles where obtained by simulation of lO'' 7 TeV beam 
protons. 

5.3 Apertures 

The physical aperture of the real optical elements has also to be taken into account. If a particle hits a 
collimator or a wall of the primary vacuum system, it is assumed to be lost. Each time a particle enters or 
leaves an optical clement, a test is performed to determine if its position matches the element acceptance. 
This is illustrated in Fig. 7. Protons from the LHC beam 1, around the IPS have lost some energy, which 
leads to a deflection with respect to the nominal beam position in (0, 0) in the transverse plane. The 
aperture shape of the MB.B9R5.B1 optical element is the rectelliptical hatched area. The proton energy 
has been chosen to 110 GeV, in order to match the element limiting energy acceptance. Due to the initial 
smearings, some protons (in black) still pass through the optical element without being stopped. On the 
other hand, some of them hit the limiting aperture shape and are tagged as stopped protons (in red). In 
such a case, the particles are considered as lost from the beam, and not propagated any further. 
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Figure 4: Top and side views of both LHC beams around the IP5 {first two graphs) and at the IPl {last 
two graphs). The IP is located at s = m. Beam 2 (red) flows from the right to the left and is seen before 
its passage at the IP. In turn, beam 1 (black) is seen downstream, passing from the left to the right after 
the crossing. In this graph, the bending effect of the sector magnets (S-Dipoles) has been switched off, 
thus rectifying the beam path to a straight line after 250 m. The difference between the IPl and IPS is 
clearly seen, as the crossing planes are perpendicular to each other. Note: In the LHC absolute frame, 
the position of the IPl is at s = m and the IP5 is at s = 13329 m. 
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Beam 1 profile at 220m 
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Figure 5: LHC beam profiles in the transverse plane at 220 m and 420 m from the IPS, for both beam 
1 (above) and 2 (below). The graphs in the left column plot the horizontal and vertical position profiles. 
The graphs in the right column show the horizontal and vertical angular profiles. These plots are obtained 
by propagating ten thousand 7 TeV protons through the LHC beamlines, in the forward direction, from 
the IPS. The mean horizontal position x = —97 mm matches the half of the separation distance between 
both beams. The red ellipse delimits the 3a contour for each plots. 
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Figure 6: LHC beam profiles in the transverse phase space at 220 m and 420 m from the IP5, for both 
beam 1 [above] and 2 (below). These graphs depict the horizontal {x^Ox) and vertical (y^Oy) distributions 
obtained by propagating 10^ 7 TeV protons through the LHC beamlines, in the forward direction from 
the IP5. 
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Aperture effect of "MB.B9R5.B1" on 110 GeV energy loss protons 




Figure 7: Example of aperture check for the MB.B9R5.B1 dipole (s = 338 m), drawn in the {x,y) 
transverse plane, at the exit of this optical element. A set of protons, with a mean energy loss around 
110 GeV, has been propagated from the IP5 to the MB.B9R5.B1 dipole using the LHC beam 1 optics. 
The beam 1 is centered in (0, 0) but is not shown. The closed blue area corresponds to the rectellipse 
aperture shape of the optical element. Protons that passed through the dipole are the black ones, while 
protons hitting the walls are tagged in red. The first graph is scaled for the full aperture shape, the second 
one has been zoomed in around the region of interest. 
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6 Very Forward Detectors 



The physics program at the LHC has been extended to processes in which isolated particles at very 
large (pseudo-)rapidities are produced [6] in difFractive and photon induced interactions [7]. To study 
those processes, several very forward detectors arc being prepared for installation along the bcamline. 
This section will discuss in particular two scenarios, differing by their distance from the IP (s) and their 
minimal horizontal distance (x) from the nominal beam position. The assumed location of the first detector 
is (s =220-224 m, x = 2000 yum), and of the second one at (s =420-428 m, x = 4000 ^m). No hypothesis 
is made on their detection efficiency or their resolution, unless quoted. We consider here VFDs providing 
2D-measurement {x and y coordinates), each consisting in fact of two stations separated by 4 and 8 m 
as a lever arm for the angle measurement, with no magnetic element in between. In the following several 
examples of studies performed using Hector are given. 

6.1 Acceptance 

Using Hector's aperture description, it is possible to identify the characteristics of the protons that will 
hit the VFD. The exchange of a photon or a pomeron, leaving the proton intact, results in a proton energy 
loss (Eioss) and a scattering angle, directly linked to the four- momentum transfer squared (t). In case of 
the photon-induced interactions —t is equal to the virtuality of the exchanged photon. The acceptance 
windows of the VFDs can be computed by performing scans in (Eioss, t) and computing the probabilities 
of reaching the detectors. The Fig. 8 shows the contour plots of the detectors acceptance, in this {Eioss, t) 
plane. It is obvious from these graphs that the VFD acceptances mostly depend on Eioss, and have a 
very small sensitivity in t, within a large t range. For VFD at 220 m, an area forbidden by kinematics is 
visible, for low |t| and high Eioss- Corresponding profiles at fixed virtuality are shown at 220 m and 420 
m (Fig. 9). 

6.2 Chromaticity grids 

Once the acceptance windows of very forward tracking detectors are defined, it is interesting to see match- 
ing between the proton variables at the IP and those measured by VFDs. Depending on their energy 
and angle at the IP, forward protons will hit the VFDs at various positions. Drawing iso-energy and 
iso- angle curves for a set of sample protons produces a grid in the measurement related variables, (xi, X2) 
or {xi, 9). Due to optics of the LHC beamlines, the grid unfolds itself in a much clearer way in the latter 
plane, and is almost invisible in the former one. The energy dependence of the transfer matrices implies a 
deformation of the grid - without such a dependence, the grid would be a parallelogram. One should note, 
that uncertainty of the transverse position of the proton vertex at the IP results effectively in smearing 
the chromaticity grids. Anyway, these chromaticity grids provide a straightforward tool for unfolding the 
energy and angle at the IP of the measured particle. The grids in Fig 10 were calculated in the energy 
range accessible to the VFDs. Finally, the correlation between the vertical angle at the IP and the vertical 
coordinate measurement is shown in Fig 11. 

6.3 Irradiation levels 

The total diffractive cross-section at the LHC is very large, resulting in a high rate of diffractive protons 
hitting the VFDs. As a result, it causes extremely high irradiation levels, requiring in turn usage of very 
radiation hard detectors. As an illustration, the number of hits and their respective map, due to the 
pp pX processes has been investigated (Fig. 12). In this approach, we neglect other radiation sources, 
like beam halo or secondary particles coming from interactions in the beam pipe for instance. 

The pp pX processes are generated by Pythia [8] 6.2.10 (single diffraction, process number 93). 
The number of hits are normalized by year by square centimeter, assuming an integrated luminosity of 
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L = 20 fb"^ From the 10^ generated events, 17 % have been detected at 220 m, and about 12 % at 
420 m. This simulation shows that the hottest spots in the detectors are locahzed with areas of only a 
few square mm, reaching proton annual proton fluencies of more than 4 x 10^* protons/cm^. 




Figure 8: VFD proton acceptance for the LHC beam 1 (left) and 2 (right) around the IPS. These VFDs are 
located at (s = 220 m, a; = 2000 fim horizontal position ; above) and (s = 420 m, a; = 4000 fim horizontal 
position ; below). This map shows contours of 25%, 50%, 75% and (plain curve) 100% acceptance. The 
acceptance is roughly rectangular, i.e. independent of t. At 220 m, the missing triangle in the acceptance 
at low \t\ and high energy loss is prohibited by kinematics. This area corresponds to non-physical protons. 
The stair-stepping effect of the lower border of this triangle is only due the binning of the graph. At 
bottom right corner of each graph, the angular kick coming from the large momentum transfers leads to 
an increasing smearing of the graph lower edge. 
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Figure 9: VFD acceptance as a function of energy loss, for two fixed virtualities (Reminder: 
See previous figure for more details. 
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Chromaticity grid at 220 m \ Chromaticity grid at 420 m 
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Figure 10: Chromaticity grids: iso-cncrgy anii iso-anglc lines for the VFDs at 220 m (left) and 420 m 
(right) away from the IPS, for the LHC beam 1. The graphs show the positions xi and X2 of protons, given 
the energy loss [0; 1000] GeV and [0; 100] GeV, respectively, and the angular kick [0; 500] fiiad. The energy 
dependence of the transfer matrices induces a deformation of the grid, worsening the reconstruction power 
at higher angles. At 420 m, the grid completely squeezes to a single line. The grids unfold themselves 
better in the x - x' plane. 
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Figure 11: Shift in vertical coordinate (yi and 2/2) for VFDs at 220 (left) and 420 m (right) from the IPS, 
due to the vertical angle 6* at the IP. 
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Hits in VFD at 220m (L=20 fb"*) 




Hits In VFD at 420m (L=20 fb"*) 




Figure 12: Irradiation levels due to pp pX processes, of VFDs located at s = 220 m (above) and 
s = 420 m (below) from the IPS. The horizontal position of the detector edge is respectively x = 2 mm 
and a; = 4 mm from the nominal position of the center of the corresponding beam. The fluence is given 
per year and square centimeter. The integrated luminosity of L = 20 fb~^ was assumed. For this analysis, 
Hector has propagated protons with 4-momentum generated by Pythia 6.2.10 (process number 93). 
The irradiation levels can locally exceed 1 x 10^^ protons per square centimeter. 
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6.4 Reconstruction 



One of the main physical motivations for Hector is the reconstruction of the event kinematics. For 
instance, if a beam particle has exchanged a photon at the IP, one could reconstruct photon's energy (E) 
and virtuality (Q^). The particle energy at a given position in the beamline is obtained from the measured 
particle position and angle within the matrix formalism by solving these equations: 

( Xs = asXo + bsXQ + dsE 
\ x'g = asXo + Psx'q + -isE 

The transfer matrix of the beamline yields the coefficients a, &, etc. The introduction of an energy de- 
pendence on the strength of optical elements refines the transfer matrix, becoming a function of E: as{E), 
bs{E),... This dependence will introduce non linearities. Let us consider two reconstruction methods for 
the energy. In the first, so-called, trivial method only the dispersion D ^ dg is considered. 

The calculation of the particle energy using a beam transfer matrix corresponding to the nominal 7 
TeV energy (i.e. assuming constant D) obviously induces some errors. Therefore, in order to compensate 
for the energy dependence of D, a polynomial fit to the average (horizontal) position at a given detector 
location was done as a function of the particle energy. Then, for a fixed four-momentum transfer (or 
virtuality Q^), the energy resolution SE is almost independent of energy itself (see Fig. 13), apart from 
the large values for the detector at 220 m. At the zero-degree angle the energy resolution is about 3 
and 1 GeV at 220 m and 420 m, respectively, but the higher is the scattering angle (or the virtuality) , the 
worse is the energy resolution, as shows Fig. 13. So, if no correction is made due to non-zero scattering 
angle, the average energy resolution for a typical of 0.3 GeV^ is worse, about 5 and 2 GeV at 220 m 
and 420 m, respectively. The trivial method is not much sensitive on the position resolution of VFDs - 
only if resolutions are close to the horizontal beam size at a given location (70 and 250 /xm, respectively), 
they would change significantly the energy resolution. However, in this method neglecting the terms UgXo 
and bsXQ leads to significant sensitivity to the non-nominal values of the average vertex position and beam 
direction (tilt) at the IP, as shown in Fig. 14. While the average vertex position can be very well measured 
using the central detectors, the beam tilt is more difficult to control - it should be known to better than 
10-20 /irad to avoid causing a significant bias. 

Once the energy has been reconstructed (Fig. 13 and 14), the scattering angles at the IP can subse- 
quently be computed, after applying all the usual smearings on other variables. As it can be seen in Fig. 
11, the reconstruction of vertical angle is really good both at 220 and 420 m. Even for very modest vertical 
position resolutions of 50 /im, the reconstruction resolutions are already limited only by the vertical beam 
divergence of about 30 firad. This reflects the fact that all dipole fields are horizontal around the IP5, 
and it is effectively a simple 'geometrical' measurement. Unfortunately, as expected, the situation is much 
worse in the horizontal plane, where the errors and biases are large and increase with energy loss (see 
Fig 15). However, in the case of the 220 m detectors it seems possible to correct for the bias and recover 
the angular resolution close to the ultimate 30 /irad. It is much more difficult at 420 m from the IP - 
this means that 6^ at the IP cannot be well reconstructed at 420 m with the trivial energy reconstruction 
method, and possibly could be only used to apply selection cuts. 

In the second (advanced) method one reconstructs the energy and scattering angles using position 
measurements at two detector stations at the same time. This requires to solve the equations for Xs at 
both detectors for E. The Qg, hg and dg coefficients of the transfer matrix depend on energy with rather 
complicated shapes, as they are products of many magnets matrices. One efficient way to get those is then 
simply to fit each coefficient as a function of energy. Various fitting functions were tried, but a quadratic 
fit proved to be sufficient to avoid any visible bias or resolution degradation. 

Using these fitted coefficients, one can easily get a formula for xq as a function of x at both detectors 
and of the energy. For each pair of detector x coordinates, the used method is to numerically find the 
root (a;o ~ 0) of the formula to get the energy corresponding to a xq = 0, and thus neglect the interaction 
point transverse extension. 
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This method allows to get reconstructed energy independently of the angle of the particle at the IP. 
Energy reconstruction resolutions does not degrade anymore if the particle transverse momentum pT rises, 
but the price to pay is that the detectors resolutions become critical. In particular, it means that the 
uncertainty of the reconstructed angle at a given detector location should be better than the beam angular 
divergences there (6 and 1.5 /xrad, respectively). The transverse momentum can thus also be computed 
from the matrix coefficients once the energy has been reconstructed. 

Distributions of reconstructed energy and transverse momentum for detectors at 220 m and 420 m 
are shown in Fig. 16. One can see that the reconstructed energy resolution stays very good even with 
non-negligible initial particle pt- Figures 17 and 18 show the pt and energy dependence of resolutions on 
the transverse momentum and energy, for various detectors resolutions. As expected, the energy resolution 
is independent of px but is sometimes sensitive to the energy as expected from the chromaticity grids of 
Fig. 10. The effect of the detector resolution, which was absent for the simple reconstruction method 
described previously, is now clearly visible even for an excellent resolution of 5 /im, especially for detectors 
at 220m from IP. 

In summary, the proton angular distribution at the IP affects the most the energy reconstruction. 
In contrast, the vertex lateral distribution has negligible impact. As a result, the most important beam 
parameter, which can change run-to-run, is the beam tilt at the IP. In principle, it can be indirectly 
controlled by the beam position monitors at 220 and 420 m, but more direct tilt measurements are 
favored, as by using BPMs next to the IP, or by monitoring the direction of neutral particle (photons, 
neutrons, or neutral pions) production in the zero-degree detectors (ZDDs). 

The energy resolution squared of the scattered proton aE can be then approximately decomposed into 
four terms: 

where the nominal beam energy dispersion ctq ~ 0.8 GcV, the contribution due to the vertex spread 
(jytx ~ 0.7 GeV at 420 m and 1-2 GeV at 220 m, the contribution due to the detector resolution, 
neglecting the angular effects, cr^gj is small for resolutions better than 50 /im, and the contribution due 
to the proton non-zero angle at the IP, Uang, which is very sensitive the angular reconstruction at the 
VFDs; even for a zero-degree scattering if one neglects effects due to the beam angular divergence at the 
IP, cfang ~ 1 GcV at 420 m and « 3 GeV at 220 m. 
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Energy resolution at 220 m - trivial method 
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Figure 13: Energy resolution for the trivial reconstruction, for protons measured at 220 m and 420 m 
from the IPS. The protons have exchanged a particle taking away some energy, at a given virtuality. The 
resolution is shown as a function of the energy loss and the virtuality. The effects of some error on position 
measurement, due for instance to the spatial resolution of the detector stations (5 /im and 30 /im), are 
taken into account. The beam energy dispersion is not included. 
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Figure 14: Bias studies for the trivial energy reconstruction for protons measured at 220 m (above) and 
420 m (below). As the energy is computed using only the dispersion term of the global transport matrix, 
any non nominal values of the event vertex x* and the beam tilt 9* at the IP, results in a biased energy 
reconstruction. 
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Horizontal angle reconstruction @ 220m (100 |irad) 
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Figure 15: Resolution ol the reconstruction ol the horizontal (O^.) scattering angles, for a 100 /irad scat- 
tering angle as a function of proton energy loss, for VFDs at 220 {upper plot) and 420 m (lower plot) from 
IP5. The trivial energy reconstruction method is used, and the beam divergence at the IP is not taken 
into account. 
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Reconstructed energy (50 GeV) 
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Figure 16: The so-called 'advanced method 'allows to reconstruct at the same time the energy and trans- 
verse momentum (or scattering angle). Here examples are given for detectors at 420 m (upper plot) and 
220 ni (lower plot) from the IP, corresponding to energy loss of respectively 50 and 400 GeV, and transverse 
momentum of 0.5 GeV (uniformly distributed in azimuthal angle). Ideal detector is assumed. 
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Figure 17: Resolution of the reconstruction of the particle transverse momentum p^, as a function of 
the energy loss (left) and of the transverse momentum (right), for VFDs at 420 {upper plot) and 220 m 
{lower plot) from IP5. The advanced method, using full detector information, is used. Dots correspond to 
different scenarios of detectors resolutions, namely perfect detectors {circles), 5 /xm {squares) and 30 /im 
{triangles) spatial resolution. 
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Figure 18: Resolution of the reconstruction of the particle energy loss E, as a function of the energy loss 
{left) and of the transverse momentum {right), for VFDs at 420 {upper plot) and 220 m {lower plot) from 
IPS. The advanced method, using full detector information, is used. Dots correspond to different scenarios 
of detector resolutions, namely perfect detectors {circles), 5 ^m {squares) and 30 /im {triangles) spatial 
resolution. 
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6.5 Calibration 



Usually, the calibration of the reconstructed variables at forward detectors can be well maintained using the 
physics processes and the central detectors. At HERA, for example, the elastic p meson photoproduction 
was used where the momentum of the scattered proton could be deduced from two decay charged pions 
using the central tracking. At the LHC the two-photon exclusive production of dimuon pairs seems a good 
cahbration process, at least for the detectors at 420 m. The visible cross section is large (crpp < 7 pb), 
including the acceptance of central detectors. This should even allow for a run-by-run calibration of the 
scattered proton energy scale within a full acceptance range. Finally, the expected reconstruction power 
of central detectors is excellent for such di-muon events, giving the proton energy uncertainty about 10~^ 
per event [6]! One should note, however, that using this process it is not possible to check the angular 
reconstruction, and that it has much more limited statistics within acceptance of the detectors at 220 
m. The high-energy bremsstrahlung pp ppj is a possible candidate for the energy scale calibration of 
the VFDs at 220 m. It has relatively large cross-section, about 10 nb for > 100 GeV, and photons 
are emitted into a very forward cone, so could be detected in the zero-degree detectors (ZDDs) at about 
100 m from the IP, proposed for measurements of neutral particles produced at very small angles. A 
simultaneous measurement of the scattered proton in the VFD and of the bremsstrahlung photon in the 
ZDD would allow then for cross-calibration of these two devices. In addition, these photons are emitted 
at very small angles (typically, at less than 150 /irad), so could be used to control the beam direction (tilt) 
at the IP. 

Finally, one should add that the angular measurements could be controlled by studying the recon- 
struction of the beam angular divergence at the IP, and comparing the obtained results with the direct 
measurements at the ZDDs, or with estimates obtained using the LHC instrumentation, or derived from 
the LHC luminosity. 

6.6 Misalignments 

The misalignment of LHC optical elements could have a significant impact on the measurements with 
very forward detectors. As the deflection of the particle paths depends on their positions in quadrupoles, 
a misplacement of these optical elements implies a change in the nominal beam position. In turn, as the 
accurate position measurement with the forward tracking detectors (as well as the information inferred 
from the segmentation of forward calorimeters) is referred to the ideal beam location, changing this 
reference results in a biased reconstruction of the measured particles. 

Fig. 19 shows the impact of possible shifts (0.5 mm) and tilts (0.1 mrad) of the beamline quadrupoles 
on the energy reconstruction with VFDs at 420 m. The reconstruction assumes ideal beamline in which 
only one quadrupole at a time is separately moved. Effects even higher than 10 % could be expected. The 
closer to the IP (i.e. the further from the VFDs), the more insensitive the correction. One can partially 
correct for these effects using information from the beam position monitors, but better results are obtained 
using a physics calibration process like the two-photon muon pair exclusive production, at least in case of 
the VFDs at 420 m. 

These misalignment effects and the corrections are illustrated (Fig. 20) by the study of two-photon 
exclusive production of the SM Higgs boson (pp — > pp(77)H) with Mjj = 115 GeV. The measurement 
of the energy of two scattered protons straightforwardly yields the boson mass, by means of the so-called 
missing mass method. As a consequence, an uncorrected measurement with misalignment leads to a 
bad mass calculation. A quadrupole (MQM9R5, s = 347 m) close to the detector has been shifted by 
100 /.im, and another quadrupole (MQXA1R5, s = 29 m) close to the IP by 500 /.im; in the latter case, 
the misalignment-induced change in the VFD acceptance is visible. The limitations of the beam-position- 
based corrections are clearly visible, even assuming no systematic errors, while the muon-calibration stays 
unbiased (though only a relatively small sample of 700 dimuon events was used to get the correction 
factors). 
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Reconstructed energy for quadrupole displacement of 500|ini (E =100 GeV) 
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Figure 19: The misalignment of the LHC quadrupoles bias the energy reconstruction. The graphs show 
the bias for the reconstruction (with the trivial method) of a 100 GeV energy loss assuming the misaligned 
quadrupoles at various positions. Each element is separately shifted (500 ^m) or tilted(100 /irad), assuming 
a perfect alignment for the rest of the beamline. The impact of the misalignment can be important. Even 
a perfect knowledge of the actual beam position at the VFD (right) does not compensate for this bias, 
depending on the position of the misplaced quadrupole. 



7 Summary and outlook 

A new tool for the transport of particles in beamlines was presented, exhibiting high performances thanks 
to the transfer matrix model. The tool was validated with the Mad-X software, on the LHC beamlines. 
Much information is produced from the raw transport of particles and study of the transport up to forward 
detectors is subsequently possible. The reconstruction of variables at IP is tested as well as the impact of 
the misalignment of beamline optical elements on this reconstruction. 

Hector is being ported in the software of the CMS experiment. Finally, further studies on the 
reconstruction methods, effects of the beamline misalignments and the related corrections will be led. 
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Misalignment impact on hHiggs mass reconstruction 
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Figure 20: Illustration of the cfFects in the energy reconstruction due to the misalignment of LHC 
quadrupoles. The graphs show the reconstructed Higgs boson mass in the two-photon exclusive pro- 
duction, using energy of two forward scattered protons. In the upper plot, a quadrupole (MQM9R5, 
s = 347 m) close to the detector has been shifted by 100 /um. Misaligning an optical element (MQXA1R5, 
s = 29 m) close to the IP leads to a loss of acceptance (lower plot). The reconstructed values including 
the correction due to the dimuon calibration is also plotted. In brackets, the average reconstructed mass 
and its resolution are given, without including the beam energy dispersion. 
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Appendix 



A Transfer matrices 

We define the transport matrices of optical elements as the following ones. The length of the optical 
element is I and k represents its magnetic field. We added the 6th line and row in order to take the kickers 
into account. The energy dependence of the magnetic deflection of the particle is set by rescaling the k 



factor by 



.Note : Our convention is transposed compared to [5] : xq.M ~ xi instead of a-i = M.xq 



E-AE- 

SO the matrices should be transposed. The matrices should have the following units : 



\ 1 1 1 1 11/ 

For a vertically focusing quadrupolc {H_VerticalQuadrupole) , where uj{k) = l^/k 
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For a horizontally focusing quadrupole [H_HorizontalQuadrupole) , where uj(k) ~ l\/k 
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For a rectangle dipole (H_RectangularDipole) , where r{k) 
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(6) 



(7) 



Note : r * (1 — cos{l/r)) causes a numerical sensitivity. We use 2sin^(a;/2) = 1 — cos(a;) instead. The 
rectangular shape implies an edge focusing of the dipole. In the beam frame, wc apply the edge focusing 
as described in [5] : 

l^rdip— edge— focusing l^edgel^rdipl^e 
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where ^' — l/2r and 
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For a sector dipolc {H_SectorDipole), either with their bending effect on or off, where r(fc) 
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The bending can be switched off when working in the detector frame. Obviously, this should not be 
used for non nominal energy beam particles. 
For a drift {H_Drift) : 
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For a horizontal kicker (H_HorizontalKicker) , where k is the angular kick in rad 
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For a vertical kicker [H_VerticalKicker) , where k is the angular kick in rad 
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B User guide 

B.l The usual "getting started" chapter 

In order to run HECTOR on your system, first download its sources and compile the libHector library'^ ; 

meSmylaptop : ~$ wget http : //www . f ynu.ucl . ac . be/per so/xrouby/files/hector_l_3_2 . tgz 
meSmylaptop : ~$ tar jxvf hector_l_3_2 .tgz 
meOmylaptop : ~$ cd Hector 

Hector has been tested with Root versions from 4.xx. When launched from Hector's main di- 
rectory, Root will automatically load libHector (see rootlogon.C) and display a welcome prompt. In 
case you are not using Pythia, one can simply ignore it by commenting fldef ine _include_pythia_ in 
./include/H_Parameters.h*. 

meSmylaptop : ~/Hector$ make 
meSmylaptop : ~/Hector$ root 
Ready for Hector — enjoy ! 
root [0] 

As a general comment, don't forget that every class you use in your routine should be included in the 
routine file, for instance : 

#include "H_BeamLirLe .h" 
#include "H_BeamParticle .h" 

B.2 The everyday hfe with Hector 

In this chapter, we will tell you how to use Hector to perform basics tasks such as using predefined 
beamlines and designing yours, propagating particles along those, and computing some beam properties. 
Generally speaking, the detailed constructors and methods arguments for Hector classes will not be 
detailed here, as it can be easily found in the Internet reference manual^. 

B.2.1 Genesis 1 : the beamline 

The first step to use Hector is to define the list of optical elements that the particles will cross while in 
the beam pipe. This can be done in two ways : 

• Using (existing) magnets tables 

• Building the line yourself 

Tables of optical elements 

This way is the simplest. You only need a text file containing columns with the magnet name, position, 
length, strength. Apertures can also be specified in this file. The order of the columns is meaningless but 
you should use column headers with the following "codenames" : 

'^Thc Makefile is compatible for a compilation on Linux and Windows/Cygwin systems. For such Windows, just modify 
the library extension name to LEXT=dll in the Makefile. 

*If your Root version does not contain the Pythia libraries, the make command will return an error. In that case, you 
will have to remove by hand the references to libPythia and libEGPythia in the Makefile, then run "make clean" and "make". 

^This online manual, as well as all information on Hector, can be found here : http://www.fynu.ucl.ac.be/hector.html 
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• NAME : name 

• KEYWORD : type 

• S : position along the line 

• L : length 

• KOL : dipolc strength (horizontal) 

• KIL : quadrupole strength 

• HKICK : horizontal kick angle 

• VKICK : vertical kick angle 

• APERTYPE : aperture type 

• APER_X : aperture size (X runs from 1 to 4) 

The "KEYWORD" code defining the element type can take the following values : 

• "DRIFT" : no-field zone^ 

• "QUADRUPOLE" : regular quadrupole 

• "RBEND" : rectangular dipole 

• "SBEND" : sector dipole 

• "VKICKER" : vertical kicker 

• "HKICKER" : horizontal kicker 

• "MARKER" : dummy element^ 

• "RCOLLIMATOR" : rectangular collimator 
The "APERTYPE" can be : 

• "NONE" : no aperture limitation 

• "CIRCLE" : circular aperture, the radius is given by APER_1. 

• "RECTANGLE" : rectangular aperture. The x and y sizes are given by APER_1 and APER_2. 

• "ELLIPSE" : elliptic aperture. Same size definition as the rectangular aperture. 

• "RECTELLIPSE" : intersection between a rectangle and an ellipse. The APER_X parameters give 
the rectangle then the ellipse size. 

When your table is properly written, the method to create your beamline is the following one^ : 

®The drift spaces lines in the tables are not read by Hector, because it automatically fills the gaps between other elements 
with drifts of the right length. 

'^This can be used to tag special places such as Interaction points. 

*As all the methods/functions using a H_BeamLine as argument require it to have the "pointer" type, it is heavily 
recommended to declare it with that type in all the cases. 
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float length = 500; 
int direction = 1; 

H_BeELmLine* mybeamline = new H_BeamLine (direction, length) ; 
mybeELmline->f ill ( "table . txt " , 1 , "starting point") 

The two arguments of the H_BeamLine creator arc the hne length and the direction in which the table 
will be read (forwards (1) or backwards (-1) ). The fill Method requires the name of the table file, the 
propagation direction of the particles from the starting point and the name of the starting point marker. 

The meaning of these two direction parameters can seem unclear, so let's suppose you use the following 
table : 



NAME 


KEYWORD 


S 


" element 1" 


"rbend" 





"element2" 


"vquadrupole" 


10 


"starting point" 


"marker" 


20 


"elements" 


"vquadrupole" 


30 


"element4" 


"rbend" 


40 



If the H_BeamLine creator argument is 1, we read only the last three elements. Has it been -1, we 
would have read the first three. The fill argument being 1, the particles will go from top to bottom. If 
it is -1, it goes from bottom to top^. 

Examples of elements tables - taken from the LHC IPS region - can be found in the "data" directory 
of Hector as a working example. 

Element-by-element beam building 

Another way to build the beamline is to add each element separately in the beamline. This is done by 
this kind of code : 

float length = 500; 
int direction = 1; 

H_BeamLine* mybeamline = new H_BeEimLine (direction, length) ; 
float position = 10; 
float strength = -0.001; 
float length = 1; 

H_VerticalQuadrupole hvq("myquadrupole" .position, strength, length) ; 
mybeamline->add(hvq) ; 

This way is not recommended, as the other one is more convenient in most of the cases. But it still 
can be used to add markers or user-defined elements useful only in specific routines, or to design a new 
beamline automatically from Monte Carlo techniques. 

It is important to note that the H_Dipole and HjQuadrupole classes are purely abstract and should 
then never be used directly. 

Alignment effects 

After building the beamline, one can move elements around their nominal position using the following 
method : 



^It can be useful to notice that the order of the Unes does not matter as the ordering is done from the S parameter. 
However, the elements have to be on their real position (before or after) with respect to the starting point. 
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string name = "MCBXA.1R5" 
double delta_x = 50; 
double delta_y = 30; 

mybeELinline->alignElement (nernie ,delta_x,delta_y) ; 

This will displace the element called MCBXA.1R5^'^ by the specified amounts in the x and y directions. 
Corresponding angles can be changed using the similar method : 

mybeamline->alignElement (name ,delta_theta_x,delta_theta_y) ; 

B.2.2 Genesis 2 : The particles 

Using the lines : 

double mass = MP; 
double charge = 1 ; 

H_BeamParticle myparticle (mass , charge); 

one particle with given mass and charge^^ is created at the "starting point" (IP) with energy defined 
by the BE (beam energy) given in the H_Parameters.h file. The other main variables in this file are : 



• 


SBE 


: beam energy dispersion 


• 


PX : 


X position at IP 


• 


PY : 


y position at IP 


• 


TX : 


X angle at IP 


• 


TY : 


y angle at IP 


• 


SX : 


X dispersion of the beam at IP 


• 


SY : 


y dispersion of the beam at IP 



• STX : X angle dispersion of the beam at IP 

• STY : y angle dispersion of the beam at IP 

The X and y positions are always given in /xm and the corresponding angles are in /xrad. The s 
coordinate is conveniently expressed in meters. 

The particle can then be smeared^^ using the SBE, SX, SY, STX, STY parameters by the following 
methods : 

myparticle . smearEO ; 
myparticle . smearPos () ; 
myparticle . smearAngO ; 

Other useful methods for the particles at the starting point include : 

H_BeamParticle : : emitGamma(f loat Energy, float virtuality) ; 

^''Please note that in LHC tables the names include the double quotes, meaning that the name should be something like 
"\"MCBXA.1R5\"". 

'^^If not specified, the default particle used is a proton. 
^^Thesc methods use a simple Gaussian smearing. 
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which allows to simulate the emission of a virtual particle with given properties, causing an energy loss 
and an angle displacement. Another way to simulate an energy change is to use the 

H_BeamParticle : :setE(float Energy) ; 

method which only changes the energy parameter of the particle without any other effects. The effect 
of the energy loss is obviously simulated by Hector 

B.2.3 Exodus : the particle propagation 

Once the bcamlinc has been set and the particle has all the desired properties, the latter can be propagated 
along the former using the following : 

myparticle . computePath(mybeamline) ; 

This computes the positions and angles of the particle at the entrance of each element and allows to 
interpolate its position in all the drift spaces. This method is the most time-consuming in Hector and 
thus defines the speed of the program. The typical timescale is less than 4 fis per particle and per element 
(including drifts) 

Aperture check 

If one wants to check if the particle has been stopped by any aperture, the method following 

bool stopped = myparticle . stopped(mybeEiinline) ; 

should be used. It checks all the element apertures and stores the closest intercepting element from 
the IP. This element can then be obtained using : 

const H_OpticalElement* myopticalelement = myparticle . getStoppingElement () ; 
cout<<myopticalelement->getName «endl ; 

for instance to show the name of the element. 
Getting the particle information 

Once the path has been computed, it is quite legitimate to be interested on the particle transverse posi- 
tion^^ and angle at some position in s. This is achieved the following methods : 

float position = 100; 
myparticle .propagate (position) ; 
float X = myparticle . getXO ; 
float y = myparticle . getYO ; 
float theta_x = myparticle . getTXO ; 
float theta_y = myparticle . getTYO ; 

The particle can then be "propagated" to any other other position and the variables obtained at that 
new positions by the same methods. The complete path can be obtained via the ROOT TGraph class 
using the following : 

^^The energies and virtualities are always given in GeV and GeV'^. The virtuality is by definition a negative number. 
^*Tested on a 1.7 GHz Ccntrino processor running Linux. 

^^The X and y coordinates are respectively the transverse horizontal and vertical displacement of the particle. 
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TGrapli* ppatli_x; 
TGraph* ppatli_y; 

ppath_x = pi . getPathCO, color_x) ; 
ppath_y = pi . getPathCl , color_y) ; 
ppath_x->Draw("AL") ; 
ppath_y->Draw("AL") ; 

A convenient way to deal with big numbers of particles is to use the H_Beam class. It builds the beam 
in the following way : 

int number_of _particles = 1000; 
H_Beam mybeam; 

mybeam. createBeamParticles (number_of _particles) ; 

where the createBeamParticles method smears all the particle variables using the pre-defined param- 
eters. Most of the single-particle methods can then be applied to the full beam, such as computePath and 
propagate. One can access one particle of the beam using the simple 

int particle_i = 10; 

H_BeamParticle* myparticle = mybeam. getBeamParticle (particle_i) ; 

for instance to get it's properties or to use methods such as emitGamma on it. Using a H_Beam is 
especially useful when the user is interested in global properties of the beam, such as its spatial extension, 
given by the (3 functions. This is achieved in a straightforward way using the method : 

float position = 100; 

float beta_x_error , beta_y_error ; 

float beta_x = mybeam. getBetaXCposition, beta_x_error) ; 

float beta_y = mybeam. getBetaY(position, beta_y_error) ; 

which returns the value of the [3 function at the given position and the error on this value. 
B.2.4 Kings : Remarks for the courageous user 

In this section we will detail some features of Hector that most people should not use, some problems 
that will appear to only few of the users, and other exceptionally boring things that could help you in 
some very seldom cases. 

Absolute frame 

First of all. Hector is designed to work in "relative frame". This means the normal transverse position 
for a particle with no angle or displacement at the IP is always 0, and then particles with nominal beam 
energy will not be deflected by the dipoles. 

However, in some case it can be interesting to check the absolute trajectory of the beam. This can be 
achieved with Hector, if and only if you can consider all your beam elements as parallel. In that case, 
it suffices to add two lines at the beginning of your routine : 

extern bool relative_energy ; 
relative_energy = false; 

You can also displace all beam elements after a given point laterally using the following : 
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float s_start = 100; 
float x_offset = -0.1: 

mybeELmliiie->of f setElements (s_start ,x_of f set) ; 

This will displace horizontally all elements from 100m onwards from the IP, by 10cm in the negative 
X direction. 

Kickers 

One exception to the "no-deflection" rule of the previous section is obviously the presence of an initial 
angle of the beam compared to its natural propagation direction, for instance a crossing angle. This kind 
of case is usually coped with using "kickers" which deflect the beam by a given angle. These kickers are 
obviously included in Hector and are the only elements that doesn't respect the "no-deflection" law. If 
one feels more comfortable to switch this effect off, it is easily be done by including : 

extern int kickers_on; 
kickers_on = 0; 

at the beginning of the routine. 
B.2.5 Revelation : compiling and running your code 

Once you have written a routine and made sure you put it in the right place - for instance the very 
convenient HECTOR/routines directory - you only need to compile it. This is most easily done by running 
Root and doing the foUowing^^ : 

meSmylaptop : ~/Hector$ root -1 
Ready for Hector — enjoy ! 
root [0] .L routines/myroutine . cpp++ 
root [1] myprogramO 

B.3 A simple example 

As a reward for your patience, here's a little working example of Hector plotting a transverse view of 
the beam at 220 m from the IP. 

// C++ #includes 
#include <iostream> 

// ROOT #includes 
#include "TH2F.h" 

// Hector #includes 
#include "H_BeamLine .h" 
#include "H_BeamParticle .h" 

using namespace std; 

void drawProf ile { 

^^Here wo suppose that Hector is already properly installed on your computer. 
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const int NParticle = 1000; 

H_BeamLine* beamline = new H_BeamLine(l ,500) ; 
beamline->f ill("data/twiss_ip5_bl_v6.5.txt" ,1, "IP5") ; 

TH2F* hp = new TH2F( "Positions" , " " , 100 , -2 . 5 ,2 . 5 , 100 ,-2 . 5, 2 . 5) ; 

for (unsigned int i=0; KNParticle ; i++) { 
H_BeamParticle pi; 
pi . smearPos ; 
pi . smear AngO ; 
pi . computePath(beainline) ; 
pi . propagate (220) ; 

hp->Fill(pl.getX()/1000. ,pl .getY()/1000 . ) ; 

} 

hp->Draw() ; 

} 
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